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Abstract. The distances that galactic cosmic ray electrons and positrons can travel are 
severely limited by energy losses to at most a few kiloparsec, thereby rendering the local 
spectrum very sensitive to the exact distribution of sources in our galactic neighbourhood. 
However, due to our ignorance of the exact source distribution, we can only predict the 
spectrum stochastically. We argue that even in the case of a large number of sources the 
central limit theorem is not applicable, but that the standard deviation for the flux from a 
random source is divergent due to a long power law tail of the probability density. Instead, 
we compute the expectation value and characterise the scatter around it by quantiles of the 
probability density using a generalised central limit theorem in a fully analytical way. The 
uncertainty band is asymmetric about the expectation value and can become quite large for 
TeV energies. In particular, the predicted local spectrum is marginally consistent with the 
measurements by Fermi-LAT and HESS even without imposing spectral breaks or cut-offs at 
source. We conclude that this uncertainty has to be properly accounted for when predicting 
electron fluxes above a few hundred GeV from astrophysical sources. 
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1 Introduction 

The interest in the propagation of galactic cosmic ray (GCR) electrons and positrons from 
astrophysical sources, like supernova remnants (SNRs), has recently been revived, mostly 
in light of contemporary, partly "anomalous" measurements [1-6]. In the context of their 
possible explanation as exotic signals from dark matter annihilation or decay, it is usually 
believed that astrophysical sources of GCRs only lead to local power law spectra, in contrast 
to strong features from exotic sources. This however relies on a simplified picture of propa- 
gation assuming not only a steady state situation but also a continuous spatial distribution 
of sources. 

On the time and distance scales of GCR propagation, however, astrophysical sources 
can be considered discrete. Diffusion renders the fiuxes originating from point-like and from 
extended sources like old SNRs indistinguishable. Furthermore, although the details of the 
injection of particles from SNRs is not fully understood (for a discussion, see ref. [7]), the 
bulk of the accelerated GCRs is expected to be released with the onset of the radiative phase, 
that is within a relatively short time. 

For nuclear GCRs, the spectrum depends rather weakly on the exact distribution of 
sources in space and time since protons and nuclei difi'use over distances of kiloparsecs before 
escaping from the cosmic ray halo, thereby averaging over the distribution of sources on 
these scales. It is therefore admissible to neglect the small scale distribution and to assume a 
steady and continuous distribution of sources. Such a distribution function can be obtained 
by generalising information from radio or x-ray surveys of source candidates [8, 9] and is 
being used in many computations of GCR fluxes including the most well-known numerical 
[10-12] and semi-analytical codes [13, 14]. 

The propagation of leptonic GCRs is, however, dramatically different as electrons and 
positrons suffer from strong energy losses due to interactions with the galactic magnetic flelds 
and interstellar radiation fields (ISRFs). Therefore, the diffusion- loss length, i.e. the distance 
electrons and positrons can travel away from the sources without loosing virtually all energy, 
is much shorter than for protons and nuclei. In addition, this distance is energy dependent 
such that at high energies, that is above 100 GeV or so, mostly young and nearby sources 
contribute and the spectrum at these energies strongly depends on the history and spatial 
distribution of sources within a kiloparsec. Therefore, it is not admissible to neglect the 
nearby small scale distribution. 
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It has been proposed [15] to model the electron-positron flux by considering a contin- 
uous distribution of sources for distances > 1 kpc or ages > 10^ yr, and to supplement it by 
the few known nearby SNRs. However, this requires that all nearby sources are known from 
surveys which seems, at least, challenging. First of all, such studies suffer from various se- 
lection effects. Radio surveys, for instance, are insensitive to SNRs of low surface brightness 
and also young but distant ones [16]. Furthermore, for individual objects it can be difficult 
to decide whether they can contribute to the local electron-positron flux. For example, it 
has been suggested that pulsars can accelerate electrons and positrons but it is not clear if 
theses high energy particles can escape efficiently from the surrounding pulsar wind nebulae 
(PWNe). Recent observations [17], however, suggest that this may be the case. In addition, 
the prediction of a flux requires accurate distance information. Considering the uncertainty 
in the distance estimates of some of the known SNRs (for example, for RX J1713. 7-3946 dis- 
tance estimates vary between 1 and 6 kpc [18, 19]) we cannot expect the distances of sources 
yet to be discovered to be much more reliable. Most importantly however, surveys using 
electromagnetic radiation will only tell us about sources on the past light cone. The propa- 
gation of charged GCRs is however dramatically different from the rectilinear propagation of 
photons such that much older sources can potentially contribute. Consequently, relying on 
surveys we might simply miss old but nearby sources which could lead to artificial features 
in the predicted electron-positron flux (for an illustration of this, see ref. [20]). 

It thus seems difficult to determine the complete real distribution of sources. However, 
as we do not expect temporal variations in the source rate to be too strong, we can generalise 
the information about sources on the past light cone to earlier times thereby building a 
model of the statistical distribution in both, position and time. The predictions possible 
using this statistical information are however only stochastic, that is one can predict the 
expectation value for the flux from a statistical ensemble of sources. For protons and nuclei, 
the fluctuations around the expectation value are rather small but for electrons and positrons, 
the flux from a particular system in the statistical ensemble, that is a particular distribution 
of sources like the one in the Galaxy, will in general differ from the expectation value. In 
this sense our ignorance of the true distribution of sources in age and distance introduces 
an uncertainty into the predicted electron-positron flux, in particular at high energies. This 
uncertainty can in principle be investigated by Monte Carlo methods [21-23], that is by 
computing the fluxes from a large number of systems, randomly drawing the source positions 
and ages from the statistical model. Here we aim at a fully analytical calculation which will 
allow us to check against possible prejudices, e.g. the gaussianity of the probability density 
for fluxes. Furthermore, this approach provides the transparency to trace the propagation of 
model parameters into the final result. 

The paper is organised as follows. In section 2 we briefly review the usual setup for 
propagation of GCR electrons and positrons in a purely diffusive model. The calculation 
for the expectation value of the flux is outlined in section 3 and the determination of the 
uncertainties is explained in section 4. We present and discuss our results in section 5 and 
conclude in section 6. 



2 Propagation Setup 

The propagation of GCR electrons and positrons is governed by the diffusive transport equa- 
tion [24], 

^.X7-{D.V)n- — {b{E)n)=q. (2.1) 
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The second and third term on the left hand side describe diffusion, with the diffusion coeffi- 
cient D, and energy losses, with the energy loss rate b{E), respectively. We take b{E) = IqE'^, 
parametrising synchrotron radiation and inverse compton scattering in the Thomson regime. 
The right hand side denotes injection of electrons and positrons from the sources. The sub- 
dominant secondary production during the propagation of nuclear GCRs is being ignored. 
As the propagation affects positrons in the same way as electrons, the stochastic effects will 
be the same for both species, and we consider only the differential density n of electrons plus 
positrons. Furthermore, as the stochasticity effects we are interested in only show up at higher 
energies, we neglect convection and reacceleration which affect the flux below ~ 10 GeV. Fi- 
nally, solar modulation is expected to modify the spectrum only marginally at the energies 
in question, so we neglect it, too. We consider the transport in a cylinder of radius Smax 
and half- height Zmax, however only enforcing the boundary condition in z-direction, that is 

n{t, S, Zmax) = 0. 

For simplicity, we assume a common spectrum for all sources^, Q{Eq), such that the 
source term factorises, q{r,t,E) = p{r,t)Q{E). The Green's function with respect to the 
source density p(r, t) reads. 
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To satisfy the above boundary condition, we apply the "mirror-charge" method of ref. [25], 
yielding. 



G{E,T,t)= g{E,rn,t) where r„ 
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It is useful to factorise the spatial dependence into a dependence on radius s and z, 
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where the sum has been expressed in terms of the Jacobi theta function, ??3, 
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(2.5) 
(2.6) 

(2.7) 



and Zcr = 4zmax/7r. As most sources as well as our position are basically in the thin galactic 
disk, z « and x{zJ^) ^ xi^ = x{0,h- 



^In fact, the source parameters, e.g. total power output, spectral indices and cut-off energies, are in 
general different for different sources. However, this only introduces an additional source of uncertainty which 
we neglect here. In this sense, the magnitude of the uncertainty due to the stochasticity of source positions 
and ages is a lower limit on the actual uncertainty. 
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3 Expectation Value of the Flux 



The flux from a source that injected a spectrum Q{E) = QQE~"'e~^^^'''^^ of electrons and 
positrons a time t ago^ at a distance s from the observer is given by the Green's function G 
(cf. eq. 2.6) of the diffusion equation, 

j^^E) = ^G=^{7rfy'e-'"/'"QoE-^l-boEty-^e-^^^^—x(^) , (3.1) 

where c/(47r) denotes the "flux factor" for relativistic particles. The flux J oi N identical 
sources at distances {sj} and times {tj} is the sum of the individual fluxes, 

J = ^ J,(i?) = ^ {Trf)-' e-''/'"QoE-^l - h^Euy-\'^^^^ —x ( ^ ) • (3.2) 
i=i i=i ^^^^^ 

We note that with the above form of the source spectrum all kinds of astrophysical sources 
can be modelled, e.g. SNRs or pulsars. 

If the central limit theorem was applicable, at a fixed energy E the fluxes J{E) for 
different realisations of N sources drawn from the same probability density, fs,t, would follow 
a normal distribution with mean //j and standard deviation ctj, 

=^N^a = ^N{G) , (3.3) 

=|^^^<? = ^^V(G2)-(G)2, (3.4) 
were (C") denotes the moments of the Green's function G = G{E, s, t), 

(G-) = I dgfG{g)g"'. (3.5) 

As a function of the random variables s and t, G itself is a random variable with the prob- 
ability density fc- In the case under consideration, s and t are assumed to be independent 
random variables with probability densities fs and ft, respectively, and thus the joint prob- 
ability density fs^t factorises, fs,tis,t) = fs{s)ft{t). 

We assemble a realistic probability density for source distances s by modelling the 
Galaxy as a logarithmic spiral [26] and weighting it in such a way that the distribution in 
galacto-centric radius agrees with what is known from radio-surveys of SNRs [8] , see figure 1 . 
The distance distribution is then obtained by transforming to a helio-centric coordinate sys- 
tem and averaging over polar angle. For details, see ref. [20] . To implement this analytically 
complicated expression for fs, we expand it as a power series in s^, see figure 2, 

oo 

fs{s) = Y,a2^s{sy. (3.6) 
i=0 

In practice, we truncate the series after 14 terms and cut off the source distribution at Smax 
to prevent the series from diverging. We emphasise that arbitrary source distributions can 
be expanded in such a power series. 



^In the absence of a reliable time-dependent model of GCR escape from a SNR during the Sedov Taylor 
phase we assume instantaneous injection. We have checked that for reasonable lifetimes 0(10'') yr this is 
accurate for all except extremely young and nearby sources which are anyway very unlikely. 
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Figure 1. The assumed distribution of SNRs in Figure 2. The probabiUty density for the dis- 
the Galaxy; the cross denotes the position of the tance of a random SNR from the Sun. 
Sun in between two spiral arms. 



Furthermore, we assume that the probability density for the source ages, e.g. the super- 
nova rate, is constant over the time scales considered, O(IOO) Myr (electrons and positrons 
of GeV energies can diffuse over 0( 100) Myr before losing their energy but are, of course, 
subject to escape losses). Therefore, we take the sources to be equally distributed up to a 
maximum time tmax = 1/ (&o^min)i set by the minimum energy -Bmin to be considered, 

p / ,\ J 1/^max for ^ t ^ ^max j /o ^\ 

^*^*)-\ otherwise. ^^■^> 
Rewriting the m-th moment of the Green's function as an integral over s and t, we have, 

(G-) = I dgfGig)g"' (3.8) 

dt / ds/,,t(s,t)G"^(s,t) (3.9) 
Jo 

= / dtm f^{s)—G^{s,t). (3.10) 
Jo Jo 2s 

With eq. 3.6, the integral in can be performed and we substitute t for the injection energy 
Eq = E/{1 — boEt). For m = 1, this gives. 



(G) = r^^f;a2.(-m)-(^+^) r^A^y 



r ( 1 + i, ) - r(i + f) 



xQoi?o"^e-^"/''^"'-xf^') . (3.11) 

^cr V ^cr / 

with r being the (incomplete) gamma function. 
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Ignoring the cut-off, E'cut — ^ oo, and expanding the term x(^^) for i'^ <C 2^ which is only 
justified for energies above ~ 50 GeV for Zmax — 4kpc and amounts to ignoring the boundary 
(condition) in the z-direction, {1/ Zci)x{^'^) — 1/\/ttP, we can get analytical estimates for 
the moments (G™), 



2 ^^^ "^ ^^....^ , ..... Vfco(l-5) 



-3m/2 



A:! 

k=0 



where 



^ r(mf5?-j^)r(2 + i-im) 

C-ife - — y , (3.13) 

.2-7 1 \ ^-mal,^rrf^'^-^ 1 1 2 



D^k=r^^mJ-^-—je-"''^-^-U^^mJ-^-—,-l-i + k+-m,mai,,j , (3.14) 

with cT^ax = ^0(1 — '^)/(4-Do)i?^^^s^ax ^^"^ Kummer's confluent hypergeometric function 
U{a, b, z). In particular, the expectation value of G is, 

(G> = J-Vfgo(^o(i-^))-r(i^)f a. ^ 



xr(i + i) 




For a homogeneous distribution of sources in a disk around the observer, oq = 2/s 
and Oj = for i > 1, 



2 

max 



4^ 4vr V4Afeo(l - S) vrs^x ^max r + 

Choosing, for example, 5 = 0.6 and 7 = 2.2 this gives a spectrum proportional to close 
to what has been measured by Fermi-LAT [5]. 

For m = 2, however, the integral in eq. 3.15 and therefore (G^) diverge because of the 
finite probability density for having an arbitrarily close source which leads to an arbitrarily 
large flux. This divergence also occurs for cosmic ray protons and nuclei as was first pointed 
out [27] in the context of an energy independent diffusion model. 

In principle, this could be cured by introducing a minimum distance Smin for the distri- 
bution fs{s) or a minimum time tmin for the distribution ft{t). In fact, it has been argued [28] 



that the absence of very young, nearby sources justifies a lower limit tmin = l/v 47rL'(£')(TsN 
in the time integration eq. 3.9, derived from the average distance to the nearest source as- 
suming a homogeneous source density rate ctsn = ^/(^•Smax^max). This statement can either 
be factual, i.e. no flux from such sources has been observed, or statistical, i.e. the contribu- 
tion of such sources to the variance is negligible. As to the former statement, the energies 
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> few TeV at which such sources would show up are stih beyond the current experimental 
reach. As to the latter one, the divergence of the variance exactly shows that such sources 
contribute significantly in a statistical model. 

Without regularisation, the standard deviation of the Green's function therefore does 
not exist. In addition, the central limit theorem cannot be applied anymore and the standard 
deviation of the distribution of fluxes does not exist either. As we will see in the next section, 
the reason that the expectation value is finite but not the variance is that the probability 
density fcig) has a broad power law tail, oc 5^"^ with a < 1. For such generalised 
central limit theorem [29] is applicable: If the probability density fcig) behaves like 
for g — )• 00, the centred and normalised sum Xjy of independent and identically distributed 
(iid) random variables Gi converges against a stable distribution S{a, /3, 1,0,1) [30]. In 
general, the distribution function for S is not known analytically but can be calculated as 
the inverse Fourier transform of its characteristic function. To determine the parameters a 
and /3, we need to find the asymptotic behaviour of the probability density faig) for large g. 



4 Uncertainty of the Flux 

The cumulative distribution function Fcig), i.e. the probability to find G < g, is given by 

FG{g)= If dtdsmUs), (4.1) 

J JVoKg 

where 

VG<g = {{t,s)\0<s< w, < t < W, G<g}. (4.2) 
With G = G{E, s,t), G < g can be transformed into a condition on s, 



> sli^ = -eHogg + flog 



see figure 3, and hence 



^max 



{7rfr'^^QoE~"'{l - &oi^t)^-'e-^/«^-^o^*)^-')l , (4.3) 



Fg= / dtftit) / dsfsis) (4.4) 

Jo Jmax[0;Sniin{E,g,t)] 







dtftit) I -^fs{s)+ dtm —fs{s), (4.5) 



where Smin > for < t < and 

■Smin < for t ^ t^,. The probability density can then be 

obtained by differentiating. 



JG = = / dtft{t) = > aai / dte (s^iJ . (4.6) 

O-y JO 9 ^Smin ^ i^max 9 ~^ JO 



We now substitute. 



A2 = ¥LJ1e^-H'' = 1 - (1 - hoEt)^-' , (4.7) 
4-Do 
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and find 



1 1 1 



X I - log 5 + loj 



(4.8) 



foo(l-'5) 



2-7 E 



^logA^ 



where A^ = 1 — (1 — h^Et^^Y ^ ■ ^^or 5 — )• c», we anticipate A^ < A^ ^ 1 and therefore neglect 
the (1 — A^) terms, 



2 /' \2^l+' 



dA^ (A2) 



logff + log^- -logA^ 



(4.9) 



where A = ((47rL»o)/(fto(l - 5)) ^^'^ QoE'^e'^/^'^'^K With the same approximation, we solve 
2 _ n w \2 



•5min = 



A 



47rZ)o 



2/3 — i-E/ -Bcut 



(4.10) 



and we finally compute 



/g ^ — ^ 5^ a2. (^-J (2 + i)-(i+^)7r-(2+*)r(l + i,{2 + i)i6 - 1) logi? 



X E- 



5+(5-l-|7){2+i)g§(2+i)^_|(2+i)E/ij^,,^_2(2+i)-l 



(4.11) 



For 5^ — 7- 00, the non-zero term with the smallest i determines the asymptotic behaviour, 
such that for ao 7^ 0, we have 



fa - r^^?i^-^-t7Qf e-|^/^-,-V3 . 

Imax OTT^Vo 2 

The asymptotic behaviour of the distribution function Fq for (7 — ?• cxd is therefore 



(4.12) 



1 - Fo{9) - wg-'/' with y^ = ~^"iE-'-'^^Qt/'e-l^/^-^, (4.13) 

and Eg = for g < 0, see figure 4. 

The generalised central limit theorem [31] then states that the centred and normalised 
sum Xj\f, 



I ( ^ \ 



(4.14) 
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Figure 3. Parameter space for the cumulative 
distribution function Fq for fixed E and g. The 
domain 'Dc<g is limited by Smax from the top and 
the maximum of and Smin(£'7 9, t) from the bot- 
tom. 



fa 




Figure 4. Probability density fa for the Green's 
function. The dashed line shows the long power 
law tail for 5 — > oo. 



weakly converges to the stable distribution 5(q, /3, 1, 0, 1) with a = 4/3 and (3 = 1. The 
normalisation constants are 



UN = NfiG and vjy 



TTW 



3/4 



2r(|)sin(f|)^ 



Ar3/4. 



(4.15) 



The expectation value /xj for the flux on Earth is therefore still c/ {A'k)N but instead 
of the standard deviation, we use quantiles Xq of the stable distribution to define uncertainty 
intervals, 



2^95% 



(4.16) 
(4.17) 



For 5(4/3,1,1,0,1) the 2.5%, 16%, 84% and 97.5% quantiles are approximately -3.6, 
—2.6, 1.0 and 8.7, respectively. The stable distribution is asymmetric with respect to its 
expectation value x = and so are the uncertainty bands with respect to /ij. In particular, 
the expectation value is larger than most of the particular realisations of the ensemble and 
in this sense not representative which testifies to the long power law tail of fc- We therefore 
define a most likely flux through the maximum Xmax — —2.0 of the stable distribution, 




3; sin (I I 




(4.18) 



The energy dependence of both, the quantiles of J and the most likely flux Jm, is 
contained in oc w'^^'^ cx i^'^')'^3(5/4g-£;/_Ecut ^ which is harder than the expectation value fij, 
so the uncertainty is growing with energy, as expected. We note the peculiar dependence of 
the uncertainty relative to the expectation value, i.e. vn/un, on the source and propagation 
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Table 1. Summary of parameters. 



Diffusion Model 



Do lO^S cm2 s 

6 0.6 

^max 3 kpc 
^0 



10-16 GeV-is-i 



from GCR nuclear secondary-to-primary 
ratios 

ISRF and B energy densities 



Source Distribution 


tmax 3 X 10« yr 
N 4.4 X 10^ 


from £^min — 1 GeV 
from supernova rate, Smax and tmax 


Source Model 


Qo 7.8 X lO^y GeV"^ 
7 2.2 

Smax 10 kpc 


1 fit to absolute e"*" + e flux 



parameters. In particular, the uncertainty increases as QqN^^^ whereas the expectation goes 
like QqN . The relative fluctuations therefore decrease with the number of sources as N^^/'^ 
which is slower than would have been the case for the normal central limit theorem where 
aJ/^XJ (X A^-V2. 

We note that although the normalisations of the quantiles of J and of the fluctuations 
estimated from the square root of the regularised variance {{5N)'^)^^'^ in Ref. [28] are different, 
their energy dependence, vi\i oc 

^_^_35/4^ is the same. This is due to the particular choice 
for the regulator tmm which corresponds to a cut-off of the probability density at Jc . It turns 
out that Jc must itself be a quantile of the probability density and for sufficiently large Jc it 
can be shown that the energy-dependence of {{5N)'^y^'^ reduces to vn. Unfortunately, the 
definition of tmin does not make manifest which quantile (((5A^)^)i/^ corresponds to. 

5 Discussion 

As an example of the analytical results, we consider the interstellar electron-positron flux 
from SNRs distributed according to eqs. 3.6 and 3.7. The parameters used for the source and 
diffusion model are summarised in table 1. Figures 5 and 6 show the calculated expectation 
value and the uncertainty bands, defined by the quantiles of the stable distribution, for 
Ecut oo. As explained above, the spectra are unreliable below ~ 10 GeV as convection 
and acceleration have been ignored. 

In figure 5 we also show the result of 50 runs of a Monte Carlo calculation of the 
electron-positron fiux with the same source and propagation model. Indeed, the fluxes from 
individual realisations of the source distribution fluctuate between the uncertainty bands. We 
emphasise that the probability density for the fluxes is not a Gaussian but a more general 
stable function. We have checked the statistical interpretation of the uncertainty bands, i.e. 
the fluxes of 68 % and 95 % of the source realisations lie within the respective bands. As 
mentioned above, due to the asymmetry of the probability density fc around its expectation 
value //G, most fluxes are smaller than fic but the ones that are larger deviate more strongly. 

What most of the fluxes have in common is a roll-over at a few TeV compared to the 
expectation value. This well-known propagation cut-off is different from a source cut-off 
and reflects the fact that the defacto youngest source has a finite age ti^.^^. Demanding 
£^ > then implies a maximum energy 1 / (bot^^^.^) , provided that the source i is close enough 



- 10 - 



- expectation value 



- - - - most likely tlux 




Energy [GeV] 



Figure 5. Interstellar flux of GCR electrons and positrons from an ensemble of source distributions, 
see eqs. 3.6 and 3.7. The solid line denotes the expectation value for the sum of fluxes from N discrete, 
transient sources without cut-off, and the dotted line shows the most likely flux, corresponding to the 
maximum of the probability density for fluxes. The coloured bands quantify the uncertainty due to 
our ignorance of the real distribution of individual sources, in the sense that the fluxes of 68 % and 
95 % of the source realisations should lie within the respective bands. The fluxes of 50 realisations of 
N individual sources from a Monte Carlo calculation are shown by the thin grey lines. 




Figure 6. Same as figure 5, but instead of the fluxes from the Monte Carlo calculation, we show 
the electron-positron flux as measured by Fermi-LAT [5] , and from the low energy (LE) [4] and high 
energy (HE) [6] analyses by HESS. The error bars on the data are statistical only, the shaded bands 
denote the systematic uncertainty and the diagonal arrows show the error due to the scale uncertainty. 
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Figure 7. Same as figure 5 but assuming a source spectrum with cut-off energy Ecut = 20 TeV. 




Figure 8. Same as figure 6, but assuming a source spectrum witli cut-ofT energy E^ut — 20 TeV. 
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to considerably contribute to the total electron-positron flux. In most cases, a few young 
sources determine the propagation cut-off. 

In figure 6 we compare our results to the electron-positron flux as measured by Fermi- 
LAT [5], and HESS [4, 6]. The calculated expectation value fits the Fermi-LAT data but 
overshoots the HESS data, in particular those from the high energy analysis. Choosing a 
softer injection, i.e. 7 > 2.2, or imposing a stronger energy dependence for the diffusion 
coefficient, i.e. 6 > 0.6, would make the spectrum softer improving the fit to the HESS 
data but making the fit to the Fermi-LAT data worse. Within the uncertainty the model 
might still reproduce the data, although only the lower edge of the 95 % uncertainty band is 
marginally consistent with the HESS data. 

The situation can be easily improved on by reinstating the intrinsic cut-off in the source 
spectrum, i^cut- For the bulk of (mature) SNRs, the cut-off energy is naturally expected to 
lie in the TeV range [32]. Figures 5 and 6 show the calculated expectation value and the 
uncertainty bands for ii^cut = 20 TeV. The cut-off improves the fit to the HESS data such 
that all data lie in the 68 % uncertainty band. 

6 Conclusion 

We have investigated the expectation value in the electron-positron flux from a statistical 
ensemble of astrophysical sources as well as the uncertainty introduced by our ignorance of the 
distances and ages of individual sources. We have considered the transport of GCR electrons 
and positrons in the usual diffusive setup with the same power law injection spectrum for all 
sources. Their average spatial distribution is modelled combining the distribution of SNRs 
in galacto-centric radius with a model of the spiral structure. For this as for any other 
source distribution, the expectation value can be calculated as a sum of power laws in energy 
where the coefficients are given by the Taylor coefficients of the source distance distribution. 
We found that the probability density fc of the Green's function G has a long power law 
tail. Consequently, the standard deviation is not defined and the probability density for 
the flux is not a Gaussian but a more general stable distribution. We have quantified the 
uncertainty around the expectation value by the quantiles of the stable distribution and have 
also calculated the most-likely fiux. We find that the level of fluctuations between different 
members of the ensemble grows with energy as expected and that most fluxes exhibit a 
propagation cut-off at a few TeV. In the case without source cut-off the uncertainty interval 
is in agreement with data from Fermi-LAT and marginally consistent with the measurements 
by HESS. This can be improved by invoking a source cut-off, e.g. at E'cut = 20 TeV. The 
analytical formulae we have provided allow to consider different types of astrophysical sources 
with arbitrary spatial distribution functions. We emphasise that the uncertainty in the flux 
is inherent to the propagation of GCR electrons and positrons and our ignorance of the 
positions and ages of individual sources. As such it needs to be considered when predicting 
their fluxes from astrophysical sources. 
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